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1 Introduction 



Supersymmetry is a promising framework for physics beyond the standard model. For this 
reason it is important to understand its nonperturbative aspects such as confinement and the 
mechanism of supersymmetry breaking. Usually lattice regularization provides a tool to study 
field theories in the nonperturbative regime. However, a technical obstacle arises in this type of 
regularization, namely it is difficult to keep supersymmetry (although some progress has been 
achieved for some specific kind of theories, for a review see [1]). In order to avoid this obstacle for 
large- supersymmetric Yang-Mills theories (SYM), we can use another regularization method 
known as Eguchi-Kawai reduction [2] . 

Another motivation to study SYM theories in the large-A'^ limit is that they are expected 
to describe the nonperturbative dynamics of string theory O HI O [6l [7] . For instance, (0 + 1)- 
dimensional maximally supersymmetric U (N) gauge theory is conjectured to be dual to type 
IIA superstring theory on a black 0-brane background [7]. This specific example has been studied 
using Monte-Carlo simulation [8l[9l[T0]; using these numerical techniques the stringy a' corrections 
can be evaluated [inj. There have been much efforts to study the + 0-dimensional theory [3] 
numerically, too. See e.g. [Ill 1121 ^^3 1 114j. Finally, large- Yang-Mills theories are interesting on 
their own because they might be solvable analytically [15] while preserving essential features of 
QCD with = 3 (for a recent review, see [16]). 

As previously mentioned, the Eguchi-Kawai equivalence [2] can be used as an alternative 
method to regularize large- A'^ SYM. The main idea of this method is that large- A^ gauge theories 
are equivalent to certain lower dimensional matrix models. Furthermore in this prescription the 
degrees of freedom of the reduced spaces are embedded in the infinitely large matrices. A UV 
regularization can be introduced by taking the size of the matrices to be large but finite. This 
regularization, differently from the lattice one, does not break supersymmetry. With such a 
motivation, in [T7] a nonperturbative formulation of the maximally supersymmetric Yang-Mills 
in four dimensions was proposed. The authors of [T7| have considered a particular solution of the 
BMN matrix model [18], namely a set of concentric fuzzy spheres, which has been argued to be 
stable due to its BPS nature. Expanding the BMN matrix model about this background, the 4d 
= 4 SYM was recovered through the Eguchi-Kawai equivalence. 

In this paper, we provide a nonperturbative formulation of 4d AA = 1 (pure) SYM in the planar 
limit by using the technique introduced in [T7j. There are two main motivations to extend the 
results presented in [17j . Firstly, 4d AA = 1 supersymmetric theories are more interesting as a 
candidate of new physics in the LHC, and it is important to consider the M = 1 (pure) SYM as 
a simplest example. 4d AA = 1 SYM is dynamically richer than 4d AA = 4 and given that there is 
no known gravity dual of 4d AA = 1 SYM providing analytical results, numerical simulations are 
a valuable tool. Even though in principle 4d AA = 1 SYM on the lattice can be studied without 
fine tuning it is computationally very demanding, and a detailed study is difficult (for recent 
numerical studies see [19j). On the contrary, our supersymmetric matrix models would require 
less resources, and allow a better numerical analysis of AA = 1 SYM. 

Secondly, several groups are seriously studying 4d AA = 1 SYM on the lattice, using con- 

^ In the context of lattice regularization, fine tuning means adding counterterms in order to restore supersym- 
metry in the continuum limit. 
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ventional computationally demanding numerical techniques. When the results of these studies 
become available they could be used to further check the validity of the Eguchi-Kawai regulariza- 
tion. After its validity has been further confirmed, the Eguchi-Kawai construction could be used 
to analyze field theories with extended supersymmetries [T7j, which cannot be studied by using 
the lattice unless we introduce fine tuning. 

In the first part of this paper we formulate the Eguchi-Kawai reduction of 4d = 1 SYM 
by using a BMN-like mass-deformed matrix quantum mechanics with four supersymmetries. In 
the regularization of 4d = 4 SYM introduced by [T7], only 16 out of 32 supercharges are kept 
unbroken and the restoration of the other 16 supersymmetries is not obvious, although supporting 
evidence has been found in [n\ I20j . In the present case, all 4 supercharges are manifestly kept 
unbroken and hence we expect that 4d AA = 1 SYM is recovered in the continuum limit. 

In the second part of this paper we consider noncommutative super Yang- Mills theories in four 
dimensions. These appear, for example, in string theory as effective theories on D-branes with flux. 
Given that noncommutative Yang-Mills theories have a D-brane origin, they can be regularized 
by using matrix models as explained in [21]. Super Yang-Mills in flat noncommutative space is 
obtained by studying the theory in a background satisfying the Heisenberg algebra [x, y] = i6, 
which cannot be realized using finite- matrices. It turns out that the Heisenberg algebra can be 
described at finite- level by considering compact fuzzy manifolds like the fuzzy sphere embedded 
in fiat space. The flat noncommutative space is then recovered as the tangent space to these 
fuzzy manifolds. One unsatisfactory property of this prescription is that transverse directions 
are necessary for embedding the compact fuzzy spaces into flat space (for example, if we embed 
5^ into there is one transverse direction). In the field theory description these directions 
turn into scalar fields, and therefore, only noncommutative gauge theories with scalars can be 
realized in this way. In the case of supersymmetric models, this implies that it is only possible 
to regularize theories with extended supersymmetries H In this article we show that by using our 
matrix model with the background proposed in [T7] in an appropriate limit, we can regularize 4d 
M = 1 noncommutative pure super Yang-Mills. Using this construction, the transverse direction 
becomes an ordinary commutative coordinate and as a consequence we recover pure M = 1 SYM 
with no additional scalars. 

This paper is organized as follows. In ^ we review the Eguchi-Kawai equivalence. In ^2.11 
we explain its deformation, namely the "quenched" Eguchi-Kawai model, which we then use to 
formulate 4d AA = 1 SYM. In ^ we provide a supersymmetric matrix quantum mechanics, and 
applying the method explained in H2.1\ we provide the Eguchi-Kawai formulation of 4d AA = 1 
SYM. In § 13.31 we prove that this regularization does not suffer from the sign problem. In ^ 
we provide a nonperturbative formulation of 4d AA = 1 SYM on noncommutative space. In 
Appendix |A] we introduce four-dimensional M = 1 SYM on the three sphere and express it in a 
form which is convenient for our purpose. 

2 The Eguchi-Kawai reduction 

In this section we review the Eguchi-Kawai equivalence [2] . The equivalence guarantees that D- 
dimensional SU{N) gauge theory and its one-point reduction are equivalent if the global (Zj^)^ 

^If we use the twisted Eguchi-Kawai model [25] , which is written in terms of unitary matrices, we do not need 
transverse directions. However it is difficult to supersymmetrize it. 
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symmetry of the latter is not broken. In the bosonic case, however, this symmetry is broken 
for D > 2 in the Od theory; it is unbroken only above a critical volume [22]. To cure this 
problem, deformations of the Od theory, the quenched \23\ 124 1 and twisted [25] Eguchi-Kawai 
models (QEK,TEK) were proposed soon after the original oneO In these models, deformations 
are introduced such that (Zjv)'^-unbroken backgrounds become stable. However, recently it was 
found that both TEK [271 [Ml EH] and QEK [3D] fail at very large-iV - their deformations cannot 
stabilize the backgrounds completely. On the other hand, by combining quenched and/or twisted 
prescriptions with supersymmetry, the background can be stabilized [271 131j . 

In ^2.1l we review the diagrammatic approach to the quenched Eguchi-Kawai model (QEK) [24j . 
First we consider the simplest case, namely the equivalence between matrix quantum mechanics 
and the zero-dimensional matrix model, and then we proceed with the Eguchi-Kawai construction 
of the field theory on |17j . 

2.1 Quenched Eguchi-Kawai model 

We consider a matrix quantum mechanics with a mass term, 

S^d = N j dtTr - \[X^,X,f + ^Xf^ , (1) 

where Xi (z = 1, 2, • • • ,d) are N x N traceless Hermitian matrices. The covariant derivative Dt 
is given by DtXi = dtXi — i[A,Xi]. At large- A^, this model can be reproduced starting from the 
zero-dimensional model 

= f • NTr [-l[Y,X,f - l[X„X,f + ^Xf^ , (2) 

where Y and Xi are N x N traceless Hermitian matrices. We embed the (regularized) translation 
generator into the matrix Y, 

Y^-s- = diag{pi,--- ,pn), Pk = ^ ( k - ^] ■ (3) 



By expanding Y around this background, Y = Y^'^'+A, the Feynman rules of the one-dimensional 
theory are reproduced at large- A^, as we will see in the following. 
The action can be rewritten as 



2tt 



^\T.\^P^-P,)iXkh,-i[A, 



Sod = ^-N{-^\ip,-p,)iXk),^-i[A,Xkhj +Tr{--[X„X,Y + —Xt] }. (4) 



4L ^, J, • 2 



We add to it the gauge-fixing and Faddeev-Popov terms ^ • NTr {^[Y^-9-,A]^ - [Y^-9- ,b][Y, c\) . 
Then, the planar diagrams are the same as in the Id theory up to a normalization factor. For 
example, a scalar two-loop planar diagram with a quartic interaction (see FigH]) is 

d{d-l) fl 2ttN\ a {A/2ttN) {A/2ttN) 



i,i,k=l 



^ J ■ 'Til (Pi ~ Pk)"^ + ~ Pk)"^ + 



Another recent proposal can be found in 
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Figure 1: Two-loop planar and nonplanar diagrams with quartic interaction vertex. 



^(^-1) 27r 2 f^^' dp dq 1 



The essence of this expression is that the adjoint action of the background matrix can be identified 
with the derivative and the matrix elements of the fluctuations can be identified with the Fourier 
modes in momentum space. The corresponding diagram in the Id theory is 

d{d-l) , f^^^ dp f^/^ dq 1 



4 y_A/2 27ry_A/2 27r(p2 + m2)(g2 + ^2)' 

where Vol is the spacetime volume, and hence by interpreting A and A/A^ to be UV and IR cutoffs, 
those diagrams agree up to the factor ^ • Vol. The other planar diagrams also correspond up to 
the same factor. 

The nonplanar diagrams do not have such a correspondence, but in an appropriate limit they 
are negligible. In the Id theory, by taking a planar limit they are suppressed by a factor A^~2_ 
In the reduced model, they are suppressed if the IR cutoff A/N goes to zero. To see this, let 
us calculate the two- loop nonplanar diagram in FigHJ for example. It reads — ^^^~P ^ , which is 
suppressed by a factor (A/A^)2 compared with planar diagrams. 

Therefore, by taking the limit 

AT^oo, A^oo, A^O (7) 
the Id model on M is reproduced from the Od model. 



2.2 Eguchi-Kawai construction of Yang-Mills on 

Next let us construct the Yang-Mills theory on three-sphere by using the Eguchi-Kawai equiv- 
alence. The essence of QEK is to find a background whose adjoint action can be identified with 
the spacetime derivative. So, the strategy is to find a set of three matrices whose adjoint action 
can be identified with the derivative on S^. Such matrices were found in [331 IIZ]- As in § [XI we 
take the radius of the sphere to be 2//i. 

We introduce matrices Lj which satisfy the commutation relation of the SU{2) generators, 

[Li, Lj] = ieijkL^. (8) 
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Since these matrices cannot be diagonalized simultaneously, we embed them in the following block 
diagonal form; 



Li 



\ 



L 



L 



[js+1/2 



v 



(9) 



where L^''^ is a {2js + 1) x {2js + 1) matrix which acts on the spin js representation. The size of 
the matrix N is 



N = ^{2js + 1). 



(10) 



We introduce a regularization by restricting the representation space to a limited number of js- 
Furthermore we take the integer s satisfying 



T T 
< s < -, 

where T is an even integer. We introduce another integer P > T/2 and take js to be 

. _P + s 
Js - • 

The large limit is taken in the following way 

P^oo, T^oo, N^oo, T/P^O. 



(11) 



(12) 



(13) 



By using these matrices we can relate a matrix model to a gauge theory on 5^ as follows. The 
action of 3d theory can be written as 



S- 



3d 



N /2\ 



+^Xf - i^xe'^'^XiXjXk + ifih'^'^X.iCjXk) 



(14) 



where the derivative Ci is defined by ()69|) . This can be reproduced from the bosonic three matrix 
model 



S='-^[ --[Xi,X,f - i^iei.kX'X^X'' + i^X^ 



(15) 



where A"! = {16tt^ / fj,^ N P) {16tt^ / fj,^ N P is a normalization factor analogous to A/27r), by 
expanding the action around a classical solution 



Xi — —fiLi, 



(16) 
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and identifying Li and xf'^^ with the counterparts in Od theory as 



L,^\U,-\ xf^^^xf^HiiU, (17) 



and replacing the trace and the integral by trace, 

(Kl'itr Tr. (18) 



2x3 



The UV and IR momentum cutoffs are given by fiT and fi, respectively, and we will take the limit 
so that 

/i^O, fiT ^oo. (19) 

We also require fi^P — > oo so that spacetime noncommutativity disappears (see 

Finally we would like to add few remarks. First, the background is a classical solution and 
hence as long as it is stable we do not need to quench it. Second, when we take the large- limit 
fixing the IR momentum cutoff n, to suppress the nonplanar diagrams it is necessary to change 
the background to —^Li 1^ and take k ^ oo limit. 



3 4d = 1 SYM from matrix quantum mechanics 

In [T7], a regularization of 4d = 4 SYM on M x has been proposed by using the BMN 
matrix model [TB]. In this section we will generalize this regularization to the case of = 1 SYM. 
We consider the 4-supercharge matrix quantum mechanics given by the dimensional reduction of 
4d = 1 SYM to one dimension. We consider its BMN-like deformation [33] in order for the 
matrix model to have the matrices ([9]) as a solution. Then, applying a similar identification to the 
one introduced in [T7] , we obtain a regularization of = 1 SYM on M x 5^ . This regularization 
keeps all 4 super symmetries unbroken. 

3.1 BMN-like matrix quantum mechanics 

We start by considering the following matrix quantum mechanics 

So = jI dtTr Q(AXO' + \[Xi,Xjf - '-^^^''Dt^ - ]^ij^,[X,,^\^ . (20) 

Here Xi {i = 1,2,3) are N x N traceless hermitian matrices, the covariant derivative Dt is defined 
hy Dt = dt — i[A, • ], are gamma matrices in four dimensions, and ipa are fermionic matrices 
with four-component Majorana spinor index a. This matrix model is obtained by dimensional 
reduction of 4d AA = 1 SYM to one dimension, and has 4 supercharges which correspond to 4d 
J\f = 1 supersymmetry. We deform it by adding BMN-like terms |34)^l . 

S = So + Sm, (21) 
^ In [34] more general kind of mass deformations to (|20p has been studied systematically. 
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where 



Sm = ^l dtTr (^l^^^j'^^i; + i^ie'i''X,X,Xk - ^Xf^ . (22) 

The additional terms contains a "mass" parameter /x. We also introduced a constant P, which 
will be fixed later. It is straightforward to see that this action is invariant under the SUSY 
transformations 

d^Xi = -ie-fiip, 

Se^P = ((AX,)7°^-^7*^[X„X,] + i/.X,6^^V'=)e. (23) 

Here e is a time-dependent parameter 

e(t) = e-"^*^°"%o, (24) 

where a is a constant which satisfies a — P = 1 , and eo is a constant Majorana spinor. Note 
that e{t) satisfies the Majorana condition. In fact different choices of a and /3 are related by a 
time-dependent field redefinition |34j . As we will see, a specific choice of a and /? is convenient to 
see the correspondence to 4d = 1 SYM manifestly. 

It turns out that this matrix model has the fuzzy sphere solution. To see this we set ip = in 
the equations of motion 

- [Xj , [X, , Xi]] + -fie^'^iX, , Xk] - ^t^X, = 0. (25) 

Hence, 

Xi = -fiLi (26) 

is a classical solution if Lj satisfies the commutation relation of the SU{2) generators ([8]). Here we 
take Li to be a matrices given in ([9]) in order to obtain the four dimensional theory. Furthermore, 
this solution is invariant under the SUSY transformation. 



3.2 Correspondence to 4d A/^ = 1 

Now, we consider the correspondence to the M = 1 SYM on M x S''^. By expanding around 
the solution of (p6|) . Xi = —fj.Li + ai, the bosonic part of the action becomes 



N f ( I 2 

•Sbosonic = y / I 2 ~ il^lLi, At]) 



-^af + ifie'^'^aiajak - ifi^e'^''ai[Lj,ak] ] . (27) 
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We can easily see that, by formally replacing 

[L,, • ] ^ At^Mx), o, ^ Xi{x), Tr ^ {2/iif j dfig Tr, A ^ X^, (28) 

the matrix model and the field theory can be identified. Similarly, we can identify the fermionic 
part of the action. The fermionic part of the matrix model can be expressed as 

^fermionic = ^ ("^ / <^^Tr (^7°A^ + ^] " i^, ^]) " /?//V^7^^V) • 

(29) 

Hence, by taking the parameters a and (3 of the matrix model as a = | and /? = — |, the fermionic 
part of the SYM and matrix model become manifestly equivalent. 

Using the replacement ()28p . we can see the correspondence of the SUSY transformations 
defined in (|73p and ()23p . The time dependence of the parameter e is also same for the SYM on 
M X 5^ (j75p and the matrix model (j24p . Furthermore when we take the continuum limit, we have 
to scale the gauge coupling constant appropriately with the UV momentum cutoff. 

Before concluding this section, an important remark is in order. In 4d = 1 SYM, there 
is a (Ziv)^-unbroken phase [33], that is volume- independent. (For related works, see [3S1 137]). 
One may think that, because of this volume-independence, 4d = 1 SYM is related to the 
dimensionally reduced model ()20p . However, the situation is not so simple. In order for the small 
volume limit of 4d = 1 to be described by ()20p . the Z]\f symmetry must be broken [38]. If the 
Zn is not broken, the derivative and the commutator [A^, ■ ] in the covariant derivative give 
contributions of the same order. As a consequence the Kaluza-Klein excited modes and the zero 
modes as well give effective masses of the same order. Therefore, even in the small volume limit, 
we cannot simply truncate the Kaluza-Klein modes. In the original Eguchi-Kawai model, this 
problem is avoided by using the unitary matrix. However it is difficult to keep supersymmetry 
unbroken with unitary variables. This is the reason why we have used the technique introduced 

3.3 Absence of the sign problem 

It is easy to see that the regularization described above does not suffer from the notorious 
sign problem after the Wick rotation. Firstly the bosonic part of the action is real. Therefore it 
is sufficient to see that the Pfaffian (or the determinant in the Weyl representation) of the Dirac 
operator is free from the sign problem. 

When /? = 0, the determinant is the same as that of the un-deformed model and in the 
non-lattice regularization method there is no sign problem |8]. The proof is a straightforward 
generalization of that of the zero-dimensional model [12] . For the proof, the Weyl representation is 
more convenient. First let us briefiy summarize the proof in Od theory. In the Weyl representation, 
the Dirac operator Maij^/skh which is defined by 

4>ajiMaij,f3ki'^f3ki = Tripr^'lAf,,^], (30) 

reads E| 

Maij,l3kl = ^af3 i^^tik^jl ~ ^l^lj^ik) , (31) 



Strictly speaking we have to project the (7(1) part because it gives zero eigenvalue of the Dirac operator, but 
we omit it here just for notational simplicity. In the Id theory at finite temperature it is not necessary. 
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where (fi = I,-- - ,4) are Hermitian matrices and ipjtp are complex matrices with two- 
component Weyl indices. can be chosen as 

r^ = a^ (^i = 1,2,3), T^ = i-l2, (32) 

where are the Pauh matrices. Let fa,ij to be an eigenvector of M of the eigenvalue A. Then, 
noticing that the adjoint operators iV^;jj,fc« = A^^kSji - A^ijSik satisfy N^.ji^ik = -N*-ij^ki ' a'^d 
C72r^f72 ^ -(r/^)*, we obtain 

<^aa'^a'i«,/3'/fc'^|'/3 = i^aij^pkl)* ■ (33) 

Therefore, 

Maij,l3kl{(y'^^'^)(3kl = {M^ji,pik'f(3lk)* 

= A*(aV^)aii, (34) 

and hence (o"^(^^)Q,jj is eigenvector of the eigenvalue A*. Note that they are linearly independent 
and hence the determinant is written as 

detM = JJlAip > 0, (35) 

i 

where i is a label for pairs of eigenvalues. It is manifestly free from the sign problem. 

Next let us consider the Id theory. First let us consider the case when /? = 0. After the 
Wick rotation, the fermionic part of the action is 

j dtTr (^Dtij - iljaf'iX^, . (36) 

In the momentum cutoff prescription [8], we compactify the time direction with period 1/T (T is 
the temperature) and fix the gauge symmetry so that the gauge field becomes static and diagonal, 

A{t) = T • diag{ai,- • • ,a]\f), — vr < Oj < vr. (37) 

Furthermore we introduce the momentum cutoff A G Z such that 

A A-1/2 A-1/2 

X,{t) = J2 Mt) = E Mt) = Yl 4(r)e-*'^^*,(38) 

n=-A r=-A+l/2 r=-A+l/2 

where lo = 2ttT. Here n and r run integer and half-integer values, respectively. Then the fermionic 
part becomes 

■4'aji{p)Maijp,f3klq'lp(3kl{q) (39) 

where 



^aijp,f}klq 



6 

{puj - T{ai - aj)) dafsSikSjiSpq - ^ o"^^ (XfiikiP - q)Sji - X^ij{p - q)5ik^ . (40) 
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Note that M has momentum indices p and q in this case. Because yX^ij{j>)j = X^ji{—p), the 
Dirac operator M satisfies 

^aa' ^a'jii-p),l3'lk(-q)(^}' (3 = {Maijp,/3klqT , (41) 

and hence if Mip = Xip then M(o"^93^) = A*((T^99^). Therefore the determinant of M is always 
equal to or larger than zero. 

For generic values of /3, the Dirac operator is shifted by a mass term Pn ■ 1 and the eigenvalues 
are shifted simply as A + A* + = (A + Therefore the determinant is equal to or larger 

than zero also in this case. 

The same discussion is applicable in any dimension, as long as the theory is obtained from 
4d = 1, and hence with the momentum cutoff the sign does not appear. Of course in higher 
dimensions we need to use the lattice regularization and hence the positivity of the determinant is 
violated. However the sign problem is treatable in the following sense, at least in less than three- 
dimensions. Suppose that one simulate the model by using the phase quenched action Squench = 
Sbosonic ~ log I det M| , where det M is the fermion determinant. Then the effect of the phase factor 
detM/|detM| can be taken into account by the reweighting as {O) = (O • phase)g/(phase)g, 
where ( • ) and { ■ )q represent the expectation value of the original and phase quenched models, 
respectively. If fluctuation of the phase becomes large, both of the numerator and denominator in 
r.h.s. become small and the numerical error cause fatal problem. For lower dimensional theories 
which is obtained from 4d = 1 SYM, however, as one approaches to the continuum limit, the 
phase factor goes close to 1 for most of the configurations; see [9] for Id and |39^ 00] for 2d. In 
such a case, the reweighting method works. In this sense the sign problem is treatable. It would 
be nice if similar property can be seen in the three- and four-dimensional lattices. 

Although the sign problem is treatable in the sense explained above, however, at finite cutoff 
level small sign effect remains and hence in order to calculate the expectation value precisely one 
has to perform the reweighting procedure. For that, one has to calculate the fermion determinant. 
One of the advantages of the present method is that the positivity is kept exactly, the reweighting 
is not necessary and hence by using the rational Hybrid Monte-Carlo algorithm we do not need 
to calculate the determinant. This property reduces the computational cost drastically. 



4 4d = 1 noncommutative SYM 

In this section we provide a matrix model formulation of 4d AA = 1 noncommutative SYM. 
First we explain how gauge theories on noncommutative space are obtained from the large-A'^ 
matrix models [21]. Then we discuss the finite- regularization. 

Let us start with a bosonic D-matrix model 

S = -i-lY.^r[X„XA'. (42) 

The model has a classical solutior^ 

4°^ =P;. (// = !,••• ,d), =0 (/i = d+l,--- ,1)), [p^,pA=iO^,-lN, (43) 

® This solution cannot be realized at finite A'^. 
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where l/vr is an x unit matrix. By expanding the action ()42p around it we obtain the U{1) 
noncommutative Yang-Mihs theory (NCYM) on the fuzzy space M.'^ with {D — d) scalar fields. 
The construction goes as follows: let us define the "noncommutative coordinate" 
they satisfy 

[x>',x^] = -i{e-y -In (44) 

this commutation relation is the same as for the coordinates on the fuzzy space with noncom- 
mutativity parameter 6, and as a consequence the functions of x can be mapped to functions on 
the fuzzy space W^. More precisely, we have the following mapping rule: 

f{x)g{x) ^ f{x)-kg{x), 

Tr ^ 'W^d'^x, 

where -k represents the Moyal product, 

f{x)^g{x) = /(x)exp l^-'-l) ^{Q-^y'^'d }j g{x). (46) 

Using this prescription we obtain U{1) NCYM with coupling constant 



glc = 47r252/Vdit^. (47) 

Similarly, by taking the background to be A^lf^ = Ik (/^ = 1) ■ ' ' ,d) we obtain U{k) NCYM. 
The UV cutoff is A ~ ({N/k)Vd^) , and gNC should be renormahzed appropriately. 



Next let us combine the above technique with the Eguchi-Kawai prescription. We consider the 
three matrix model, D = 3, and take the background to be 

^1°^ = Pni ® Ina, -^^2°^ = Qni ^ Ina , = l^^ (g) diag{pi,- ■■ ,^^2), (48) 

where 

iPni^Qni] = -iO ■ Ini (49) 

and 

Pfc = — f - ?V (50) 



n2 V 2 

Intuitively, the fuzzy planes extending into the (xi, X2)-direction are located at each value of X3 
(Fig. [21 left). Then, in the large- iV limit (ni,n2 — > 00) we obtain two noncommutative directions 
using the previous construction and one ordinary (commutative) direction by the Eguchi-Kawai 
prescription. Note that the gauge group is U{oo). 

To realize this configuration at finite- A^, we can use the background ([9]). First let us consider 
the single fuzzy sphere of spin j, described by (2j -|- 1) x {2j + 1) matrices. By zooming in the 
north pole, i.e. by looking only close to L3 = j, we have 

[;uLi,/iL2] ~ -i/i^j. (51) 
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Figure 2: Fuzzy planes extending to (xi, X2)-direction placed along (left) can be obtained by 
zooming in the north pole of a set of concentric fuzzy spheres (right). 



Hence, we can identify fiLi and fiL2 to be p and q with the noncommutative parameter 

= fi'j. (52) 

The tangent space at the north pole, which is identified with the noncommutative plane, is located 
at 

X3 = nj. (53) 



In order to obtain the three dimensional noncommutative flat space described in ()48p we must 
consider the tangent spaces for the whole family of concentric fuzzy spheres defined in ([9]). These 
planes must be required to have the same noncommutative parameter 0, but will be placed at 
different positions X3 = ^js- In this description, the distance between the neighboring spheres /i 
turns into the infrared momentum cutoff along the commutative direction. Therefore, we have to 
take the following limit 

^l^js ^e, ^L^o. (54) 

Let us now consider how to impose this limit in the background We first define the 
infinitesimal parameter e as /.t = e. Then, jo = P will behave as 

P ~ e-^ (55) 

In order for the noncommutative parameter to coincide, the maximum and minimum of fi'^jg must 
go to same value 

f.'I^P±^)^l.'P (56) 

Hence, T must scale as 

(57) 
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with 7 > — 2. We must also send the distance /zT between the first and last tangent planes, which 
corresponds to the UV cutoff in the Eguchi-Kawai reduction, to infinity, 

A ~ /iT ^ cx). (58) 

From the previous relation we obtain 7 < — 1. In summary, by sending T and P to infinity in the 
following way: 

Ai~e, P~e~^ T'-e'^, e ^ 0, (59) 

where —2 < 7 < —1, we obtain the 3d (Euclidean) noncommutative space from three matrices. 

Similarly, we can obtain (1 + 3)-d N = 1 noncommutative super Yang-Mills by applying this 
procedure to the BMN-like matrix model. 

4.1 Stability of the background 

For the bosonic models, noncommutative backgrounds which lead to the ordinary flat non- 
commutative space are unstable [291 [271 [31] . Such an instability was originally argued in [5T] from 
the point of view of UV/IR mixing |l2]. The argument in [271 131] can be applied to the present 
case as well and hence the 4d bosonic NCYM cannot be constructed in this way, as also expected 
from the NCYM calculation [^if"! . 

In supersymmetric models this type of instability does not seem to exist. Given that the 
concentric fuzzy sphere background is supersymmetric, we can expect it to be stable in the limit 
discussed above at least at zero-temperature. It would be interesting to study the stability ex- 
plicitly by using the Monte-Carlo simulatiorP^. Given that noncommutative spaces in bosonic 
models are unstable, as mentioned previously, we expect the existence of a critical temperature 
above which the noncommutative space destabilizes. 

5 Conclusion and Discussions 

In this article we have introduced a BMN-like supersymmetric deformation of the 4-supercharge 
matrix quantum mechanics which can be obtained from 4d = 1 SYM through dimensional 
reduction. By using it, we provided nonperturbative formulation of 4d = 1 planar SYM and 
4d = 1 noncommutative SYM. These models can be studied numerically by using the non- 
lattice simulation techniques of supersymmetric matrix quantum mechanics [Hj-^i . An important 
application which we have in mind is the analysis of the finite temperature phase structure of 
M = 1 SYM. In theories that have a gravity dual, as = 4 SYM, it is possible to infer the 
existence of a deconfinement phase transition at strong coupling. This is just the usual transition 
from pure AdS to a black hole that takes place on the gravity side. Using the formulation of 
|17j , we may reproduce the transition from gauge theory. In the case of AA = 1 SYM the gravity 

Recently it has been claimed that by using Eguchi-Kawai model with double trace deformation [2^ 4d bosonic 
NCYM can be defined [43]. If so then the deformation should eliminate UV/IR mixing somehow. It would be 
interesting to study this point in detail. 

The stability of fuzzy sphere in zero-dimensional supersymmetric matrix model has been studied by using 
Monte-Carlo simulation in |44j . 

In one dimension lattice simulation also works, that is, supersymmetry is restored in the continuum limit [9]. 
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dual is not known and as a consequence it is difficult to study the phase structure of the theory 
analytically. Our nonperturbative formulation allow us to study the finite temperature theory 
numerically and to investigate the presence of a deconfinement phase transition. However a 
possible subtlety is that the background ([9]), which is necessary for the Eguchi-Kawai equivalence, 
may be destabilized before the phase transition takes place. Also, it is important to confirm that 
the Eguchi-Kawai construction works at strong coupling. For that, it is important to compare 
with the lattice calculation. 

Another interesting direction is to use the formulation to obtain insights into nonsupersym- 
metric theories. In |45j . a large- A'^ nonsupersymmetric gauge theory with a quark in the two-index 
representation has been discussed. This model is interesting because it can be regarded as a cer- 
tain large- A^ limit of one-flavor QCD, in the sense that it reduces to the ordinary one at A^ = 3. 
It was found that this theory can be embedded into 4d = 1 SYM, and a class of interesting 
quantities such as the fermion condensate take the same value as in the counterpart in SYM at 
large-A^. Our formulation would be used to obtain insights into the nonsupersymmetric theory 
through this equivalence. In this context, it is interesting to study the nonsupersymmetric theory 
itself, not necessarily with one flavor, by using the Eguchi-Kawai equivalence. For such systems, 
two kinds of the large- A^ reductions - the one on jl7j discussed in this paper and the usual one 
with unitary variables studied in [37] - are applicable. In this case a possible disadvantage of the 
former formulation is that it is necessary to put many fuzzy spheres on top of each other in order 
to stabilize the background against the repulsive force coming from many fermions. It would be 
nice if we could study the system by using the Eguchi-Kawai equivalence. 

Another interesting direction is to formulate large- A^ field theories on more complicated space- 
time using the Eguchi-Kawai equivalence. For example, from the discussion in S I2.21 one can easily 
see that, by restricting the spin of SU{2) generators in ([9]) to be integer values, one obtains super 
Yang- Mills theory on MP'^. By noticing the similarity of the ordinary QEK and Taylor's T-dual 
prescription for matrix models [l6], a generalization of Taylor's work to nontrivial manifolds [U] 
may lead to other examples of Eguchi-Kawai construction for more complicated spaces. (Note 
that the original construction [T7j has been obtained in this way). Generalization to 4d = 2 is 
also possible. For that purpose, we have to introduce a mass deformation which preserves eight 
supercharges. Such a deformation has been discussed in |34j . 

Furthermore the BMN-like matrix model has several applications. For example, it could be 
used to numerically confirm the existence of the conjectured "commuting matrix phase" in the 
BMN matrix model proposed in |48[ I49j . (Originally such a phase has been conjectured in 4d 
AA = 4 SYM on R x 5'^ [SD])- Confirming the existence of this phase is important because it would 
enable us to study AA = 4 SYM in the strong coupling regime. This model can also be used to 
study the stability of the fuzzy sphere solution in a supersymmetric setupf^. This will provide 
further intuition about the dynamic of the fuzzy sphere solution for the BMN matrix model. 
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A N = 1 SYM on M X 

In this section, we write down the action for J\f = 1 SYM on M x S"^ in a form convenient for 
our purpose [T7] (a detailed discussion of SYM on curved spaces can be found in ). We take 
the radius of the sphere to be 2/ ^. The action of U{N) SYM is given by 

5 = 1 d'xy^^^Tr + ^i^^-Da^P^ , (60) 

where X^d is the 't Hooft couphng constant, g^u{x) is the metric and g{x) is its determinant. The 
field strength is 

F^,u = d^A^ - dyA^ - i[A^, A^\ (61) 
and Da is the covariant derivative defined by 

= dai^ -i[Aa,i;] + ^UJabcl'"'^ (62) 

for adjoint fermions. The Greek indices fi, v refer to the Einstein frame and the Latin indices to 
the local Lorentz frame. 

The sphere part of this geometry has the group structure of SU (2). Given this group structure, 
there exists a right-invariant 1-form dgg~^ and dual Killing vectors £j, satisfying the commutation 
relation 

[Ci,Cj] = Kijk^k- (63) 

Using the coordinates {9,'ip,ip) defined by g = e-«'PO'3/2g-i6i(T2/2g-r?/>(T3/2^ vielbein can be 
expressed as 

E^ = —{— simpdO + sin 9 cos (fdip) , 
M 

E"^ = — (cos ipdO + sin sin ipdip) , 

E^ = - {dip + cos e dil)) . (64) 
A* 



In these coordinates the metric is given by 

dO'^ + sin^e dp^ + {dip + cos e dp f . (65) 

The spin connection u^abc can be read off from the Maurer-Cartan equation, 

* i 



ds^ = \ 



dE' -JjkE^ AE'' = 0, (66) 



and the Killing vectors are given by 



Wjjfc = ^^ijk- (67) 



A = —EfdM, (68) 
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where 

cos 9? 



= 


—i y— sin Lpde 


C-2 = 


—i ^COS (fdg — 


C3 = 





sm r ' 



sm 



simp 

(69) 



The Kilhng vectors represent a complete basis for the tangent space on S^. Furthermore given 
that the vielbeins are defined everywhere on 5'^, the indices i can be used as a label for the vector 
fields and l-forms0 

By using the Killing vectors £i, the bosonic part of the action can be rewritten as [T7] 



^boson = dt j d^^Tri]^{DtXi-iiC,AtY 



+^{£iXj — CjXif' — ^{CiXj — CjXi)[Xi, Xj] + -[Xi,Xj] 



-i^Xf + if.e'^'^XiXjXk - ifih'^''Xi{C,Xk)j , (70) 

where Xi is defined in such a way that the 1-form of the gauge field on take the form A = XiE"^, 
and is the volume form of the unit three-sphere. The fermionic part can be expressed as 

Sfermion = £^ (-^) (^) / / ^^^S^r (^V^T^^oV' + V^f (^/"AV' " i[Xi,i']) + ^/^V^t'^'V') • 

(71) 

The SUSY transformations of = 1 SYM on curved background are given by 

SAa = ii^-iae, (72) 
Si^ = ^F'^'jabe. (73) 

The parameter e is related to the isometry of the geometry. In the case of M x S^, e must satisfy 

V.6 = ^7a7^^'e, (74) 

which corresponds to the Killing spinor equation in supergravity. Actually, given two spinors e 
and ( satisfying (j74|) . Ct^'^ is a Killing vector. The solution to (|74p is given by 

e = e-W^^o, (75) 

where eo is a constant. 



This property is necessary in order to identify this index with the one in the matrix model |53) . 
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